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Abstract 

Landau's Fermi Liquid parameters are calculated for non-super fluid pure neutron matter in the 
presence of a strong magnetic field at zero temperature. The particle-hole interactions in the 
system, where a net magnetization may be present, are characterized by these parameters in the 
framework of a multipolar formalism. We use either zero- or finite-range effective nuclear forces to 
describe the nuclear interaction. Using the obtained Fermi Liquid parameters, the effect of a strong 
magnetic field on some bulk magnitudes such as isothermal compressibility and spin susceptibility 
is also investigated. 
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I. INTRODUCTION 



Objects of astrophysical origin such as soft-7 repeaters and anomalous X-ray pulsars 
are believed to be neutron stars with strong magnetic fields and constitute an example of 
the so called magnetars [1]. On their surface magnetic field strengths can be of the order 
Bmagnetar ~ 10^^ — 10^^ G. Accordiug to the scalar virial theorem [2] allowed field strengths 
could be, in principle, as big as -B ~ 10^^ G, however, a full detailed study of the gravitational 
stability condition of the maximum sustainable magnetic field strength in a star remains to 
be done. 

Descri]3tion of low temperature homogeneous fermion systems was developed by Landau 
in 1957 [3,] with the characterization of the low energy properties of the system in terms of 
a set of parameters, i.e., the Landau parameters, which determine the effective interaction 
of two quasiparticles at the Fermi surface. Since then, there have been many attempts to 
characterize nuclear and neutron matter as a Fermi Liquid. Even if the general framework 
of the Fermi Liquid Theory (FLT) is well defined, the microscopic derivation of the Landau 
parameters in the context of nuclear liquids is not yet fully solved and their values depend 
on many-body effects and details of the nucleon-nucleon (NN) interaction that need to be 
treated carefully. Microscopic calculations of Landau parameters for nuclear and neutron 
matter for realistic interactions have been subject of several investigations . These 

calculations have also been evaluated in the framework of effective interactions, which allow 
for a Hartree-Fock description of the system and for a clean definition of the quasiparticle 
excitations. In this context they have been calculated for symmetric and asymmetric nuclear 
matter [7|, |8| . 

From the astrophysical point of view the importance of FLT is due to the fact that it 
constitutes an excellent tool to study the properties of matter in the interior of neutron 
stars where the low temperature approximation is relevant. In the solid outer crust, the 
picture of non-superfiuid neutrons is adequate, however, at higher densities in the crust one 
should consider a finite pairing gap energy [9|. In this non-superfiuid scenario one could 
approximate, at these densities, the fermion system in the very neutron rich astrophysical 
object as formed by pure homogeneous neutron matter as that described by the FLT. Due to 
the presence of strong magnetic fields in these objects the possibility of magnetized matter 
must be taken into account. 
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In a previous contribution we analyzed the effects of a strong magnetic field of astrophys- 
ical origin in a neutron gas [l^. There, we found that using different parametrizations for 
the non-relativistic effective NN interaction such as Skyrme and Gogny forces, in the con- 
text of Hartre-Fock calculations, a net magnetization is energetically allowed. One should 
notice that effective nuclear interactions such as Skyrme forces predict that even in the ab- 
sence of a magnetic field an spontaneous magnetization at sufficiently high densities may 



arise 



Q, Q, m Q. 



However, modern calculations using realistic NN potentials such as 



the Auxiliary Field Diffusion Monte Carlo (AFDMC) metho d |l5l| . lowest order constrained 



variational (LOCV) method [l^, Brueckner-Hartre-Fock 
or relativistic Brueckner-Hartre-Fock 



17, 181, relativistic mean field 



20|, which have studied the energetics of spin po- 



larized neutron matter in the absence of magnetic fields, seem to prevent an spontaneous 
ferromagnetic phase transition and therefore there is the tendency to consider this fact as 
a pathology of the model used for the effective interaction. Notice also that the modern 

n 

finite range Gogny effective interactions do not predict this instability p^. In any case, the 
transition into a ferromagnetic state could have important consequences for the evolution of 
a protoneutron star, in particular, for the spin correlations in the medium which do strongly 



affect the neutrino opacities inside the star 
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22l |. For lower densities, as those in the 
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26| 



outer crust, microscopic calculations for the non-homogeneous pasta phases 
and the spin response in the presence of strong magnetic fields need to be fully calculated. 

Therefore it is interesting to study neutron matter in the presence of a strong magnetic 
field or with net polarization. As a first step in the description of the system, besides its 
energetics, which has already been discussed in a previous publication [10], we will study 
the Landau parameters of the polarized neutron matter, characterized for having two Fermi 
spheres corresponding to the two possible spin orientations of the neutrons. This study re- 
quires the use of an extension of the FLT to a two component system, which is presented in 
section [TTl together with the derivation of the Landau parameters for two parametrizations 
of the nuclear interaction, namely effective Skyrme and Gogny forces. Using a multipolar 
expansion of the particle-hole matrix elements we will examine the Landau's Fermi Liq- 
uid parameters up to dipolar order. Other magnitudes such as isothermal compressibility 
and spin susceptibility will be also analyzed. Results will be presented in section llllt and 
summary and conclusions will be given in section IIV[ 
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II. POLARIZED NEUTRON MATTER AS A FERMI LIQUID 



In this section we briefly review the basic formahsm of the normal FLT 27| for a ho- 
mogenous neutron system in the presence of a strong magnetic field. 

We consider a uniform magnetic field in the z-direction, B = Bh which also defines the 
quantization axis for the spin. In such a fermionic system particles can have spin projection 
on the z-axis, a, which can be either a = +loro" = — 1 for spins aligned parallel or 
antiparallel to the magnetic field, respectively. Let us remind here that for a neutron the 
magnetic moment is antiparallel to the spin. The neutron number density, p, is the sum of 
spin up (+) and down (— ) particles , 

P = P+ + P~- (1) 

These densities define the Fermi momenta, i.e., the Fermi surface of each fermion component 
with spin projection cr. At zero temperature the number density for each component is given 
by (we set ^ = c = 1), 

In the neutron system the net magnetization density is defined as, 

m = PnAp, (3) 
-1.9130427(5)pAr is the neutron magnetic moment in imits of the nuclear 



where /i„ 
magneton 



28l | and A is the spin excess or polarization of the system. 



A = P-±^. (4) 
P 

The total magnetization of a given volume is then M = J mdV. The relevant thermody- 
namical potential to study neutron matter at zero temperature under an external magnetic 



field, H, is the Helmholtz free energy, F, defined as 291] 



F = E- HM, (5) 

where E is the energy of the system. Note that previously we have used B to designate the 
magnetic field strength, however, the total magnetic field is the sum of the external magnetic 
field and the induced magnetization, B = H + 47rM. In this work we will assume that the 
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ratio I H/ B \ will always be close to unity. In what follows we will keep the notation using 
B to designate the total magnetic field strength. 

In the context of FLT it is assumed that ph excitations happen around the Fermi surface, 
so, this treatment is valid, in principle, for the low temperature regime where T << Tp^- 



27(1 . Then, there is a one-to-one correspondence of states in the interacting system with 



states in a free Fermi gas. As particles with given momentum k and spin projection a are 
added addiabatically to the system, an eigenstate of the real gas is obtained and, thus, a 
distribution of quasiparticles, nk^^- 

Since the lifetime of the quasiparticle varies inversely to the square of the departure of 
its energy from the Fermi energy 27|] we will consider that the ph perturbation mechanism 
takes place around the available polarized Fermi seas. Asuming small deviations of the 
distribution function of the a-polarized quasiparticles in the plasma, rik.o-, with respect to 
the ground state distribution, we have 5nk,cr = n\^,a ~ "^ko-- 

In the FLT the variation of the free energy, 5F , can be written up to secon order in (5nk,o- 

as, 

= Xl(^k,<7 - /inO-5)(5nk,<x + ^ fk,a,k'ySnk,aSnk'y + 0{6n^). (6) 

k,o- k,k',(T,cr' 

As seen on [10|] the change in free energy under the presence of a magnetic field is due to 
the energy of (anti)alignement of spins in the field and to the change in the quasiparticle 
distribution function. The single particle energy can be obtained as a functional derivative 

dF 



as well as the quasiparticle interaction coefficients 



d^F 

/k,cr,k',o-' = 7; • (8) 

dnk,<xC'^k',<7' 

In the general case of polarized neutron matter, the ph matrix element can be written as 



3G 



3l| 

Vph{(li, q2, q) = (q + qif^i, qicrsi ^Iq + q,2(J4, q2cr2), (9) 



where we keep the same notation as defined in l30|, IMI] to designate ki = q + qi, 

k2 = q2, ks = qi and k4 = q + q2 as the participating momenta of the ph interaction with 
initial momenta qi,q2 and q is the transferred three- momentum. Assuming that the ph 
perturbation mechanism takes place around each Fermi surface, each quasiparticle incoming 
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momentum, qi, must be replaced by the corresponding polarized Fermi momentum kp^^ and 
the momentum transfer is supposed to be small q ^ 0. 

When computing the quasiparticle interaction matrix elements we must note that they 
only depend on the Fermi momentum of each polarized component and the relative angle, 9, 
of interacting quasiparticles three- momenta qi, q2. For a non polarized system it is usually 
written as an expansion in Legendre polinomials |27| 



Vph = Y.[fi + 9i(Ti.a^]Pi{cose), (10) 



Previous works 



/=0 



31 



32| have studied a non polarized pure neutron system using a 
ph interaction at the monopolar and dipolar (/ = 0, 1) level. In that case, it is convenient 
to define dimensionless parameters Fi = J'iNq and Gi = QiNq where A^o is the quasiparticle 
level density at the Fermi surface at T = 0. For non polarized neutron matter of density 
p = gkp/{Q7r'^) is given by A^^o = ^^^^2^ where g = 2 is the spin degeneracy and m* is the 
effective mass at the Fermi surface. 

Considering the two possible spin orientations a = ±1 the polarized interaction matrix 
elements can be written using coefficients depending on the polarizations involved, //'^''^'*. 
On a non polarized system they fuUfiU the following relation, 

h = , (11) 

f((T,cr) _ n{cr,-cr) 

91 = ^ f • (12) 

In the same way, as mentioned for the non polarized case, one can define dimensionless 
coefficients in the form Fi'^''^ = //'^''^ V^AW/Vq^. where A^oo- is the quasiparticle level density 
at the Fermi surface of each component at T = 0, Nqo- = '"^ . The Landau coefficients for 
the non polarized case can be recovered from the general polarized coefficients as A — 0, 

Ft' - ^ y (5 = 0), (13) 

G|'> ^ . ^ (B = 0), (14) 



where we have defined Fjf^ and Note that some works define combinations of coeffi- 



cients in a different way [33||- Contrary to the non polarized case, now the Fi''^ interaction 



displays a 2 x 2 matrix structure depending on the spin projections considered. To size the 



importance of the polarization in neutron matter it is useful to consider a ratio of these 
coefficients {Fi'^\ G\'^^) with respect to the B = case. We define 

= \ — (15) 

and 

_ 2Gr^-Q(g = 0) 

" |GK5 = o)| • ^^^^ 

In this work we will concentrate on R~pi., Rqi which correspond to the energetically favourable 
dominant population fraction with magnetic moments (spins) aligned parallel (antiparallel) 
to the magnetic field. 

In this work we will be interested in retaining just the terms with / < 1. Notice that 
constraining results just to dipolar multipolarity is exact for the zero-range Skyrme forces 
but for the finite-range Gogny forces the coefficients at higher multipolarities are not zero 
although they decrease rather rapidly. 



Within the context of the Landau FLT, additional static physical quantities of interest 3J, 



35| can be obtained. We have derived the expressions for these quantities for the polarized 



neutron matter case. The quasiparticle effective mass is related via Galilean invariance with 



the dipolar matrix elements 271], 



K/m = 1 + ^iVo.[/^^^ + {^)ft'-% (17) 



The isothermal compressibility, K = 9^, is related to the study of the variation of the 
pressure and density profiles for a neutron star. At zero temperature can be written as, 

^ = - E 1^(1 + NoAft""^ + i^)ft'-''^])- (18) 

The thermodynamical magnitude that characterizes the change in the magnetization of the 
medium when an external magnetic field exists is the spin susceptibility, x- It is obtained 
by considering the change of the quasiparticle distribution functions for the up and down 



components 



m = xB = fin{Sp+ - Sp-). (19) 
For a polarized system we have obtained 

^ = E r^f^^ — r^- (20) 

^ l + NoAft'''^-i-§^)ft"'^]) 



The bulk magnitudes, K and x fo^^ polarized neutron matter are affected by a factor involving 
the quasiparticle interaction matrix elements and the level densities with respect to the 
polarized free Fermi gas value for each population component. 

A necessary ingredient in the calculation of the free energy functional in the neutron 
system, F, is the nuclear interaction. In the next subsections of the present paper we 
consider either zero-range Skyrme or finite-range Gogny forces as illustrative examples of 
effective nuclear interactions. Using the non-relativistic Hartree-Fock approximation, the 
free energy per particle in neutron matter under the presence of a magnetic field B can 
be calculated minimizing F as a function of B for given thermodynamical conditions to 
obtain Amj„ (polarization at the minimum). As a rule of thumb, for the maximum magnetic 
field strength considered in this work, B = 10^^ G, a A^m < 40% is allowed for densities 



0.5 < p/po < 3.3 at zero temperature [10 1. 



A. Skyrme force 

We have considered, in first place, the phenomenological Skyrme interaction that appears 
in the literature under the rather general form js?] 

F^^^^"^^(ri, r^) = to (1 + XoPn 5(r) + (1 + x.P^ [k"S{r) + 6{r)k\ 

+t2 (1 + X2Pn k' ■ 5(r)k + (1 + x^Pn p"(R)5(r), (21) 

o 

where r = ri — r2, R = (ri + r2)/2 and k = (Vi — V2)/2z is the relative momentum acting 
on the right and k' its conjugate acting on the left. P"^ is the spin exchange operator. Note 
that we have omitted the spin-orbit term not relevant for the total energy of homogeneous 
systems. Its effect on the RPA response function has been shown to be sizeable only at 



values of the momentum transfer much higher than the Fermi momentum [38 1. 

As widely known, this effective interaction allows for a good reproduction of finite nuclei 
observables and bulk matter equation of state (EOS) relevant to neutron stars js^]. In 
the Skyrme model the ph interaction matrix elements retain the density and polarization 
dependence and their multipolar decomposition gives contribution up to / = 1. We have 
derived analytical expressions for both parallel and antiparallel coefficients for an arbitrary 
polarization in neutron matter. 
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The monopolar terms can be written as, 



(22) 





- xi) + h{i + X2)] (4,. + 4,-0 , 



(23) 



and the dipolar terms can be written as, 



1 



-t2(l + a;2)4. 



-[tl(l - Xi) + t2{l + X2)]kF,a^F-a- 



(24) 
(25) 



In this paper we use the SLy7 parametrization [37[ of the Skyrme force which not only 
provides good values for binding of nuclei but also a neutron matter EOS in agreement with 
microscopic calculations obtained using realistic interaction and giving values of maximum 
neutron star masses around 1.5Mq. In Table |T] we summarize the values of some observables 
for symmetric nuclear matter (SNM): saturation density, po, binding energy for symmet- 
ric nuclear matter, a^, symmetry energy, a^, and incompressibihty modulus, Koo, for the 
effective interaction models used in this work. 



TABLE I: Values of some observables for symmetric nuclear matter in absence of magnetic field 
with the Skyrme and Gogny forces considered in this work 



33, 



Model 



SLy7 
DIP 



/Jo(/m-3) (MeV) a, (MeV) a, (MeV) 



0.158 
0.1737 



229.7 
266 



-15.89 
-16.19 



31.99 
34.09 



B. Gogny force 



The Gogny force has also been extensively used in the literature. In this paper we use the 



of pure neutron matter 



DIP parametrization 39 1 w hich allows for a good description of both finite nuclei and EOS 



13| . It includes a sum of two Gaussian shaped terms that mimic 



the finite range effects of a realistic interaction in tlie medium. Usually, it also contains a 
density-dependent zero range term. The interaction potential is written as: 

i=l 

t3.(l + a;3,P'^)p"^(5(r)}, (26) 

where P"^ and P"^ are the spin and isospin exchange operators. The values of some observ- 
ables for SNM with the Gogny DIP force appear in Table I. We have derived the expressions 
of the monopolar and dipolar Landau parameters with the Gogny interaction: 

t3i(l - X3^)[2aip'''-\p - p„) + ai{ai - l)p"'-V+P-], (27) 



= Y.{w.- H,Wp\ 

1=1,2 



3/2 

kp{a)kF{—cr) 

%(1 - X3i)[(«^ + l)p"' + aiia^ - l)p°»-V+P-], (28) 



and for the f\ terms, 



i=l,2 '^F,(T '^F,aH' 

2 



1 + Tr-2)e-''^^^'')], (29) 



(1 + _)e-i(^^-+'=^.-)''^']. (30) 

It is worth mentioning at this point that, once the ph interaction matrix elements have 
been calculated, the response functions of neutron matter can be obtained along the lines 
described in [38f|. 



III. RESULTS 



In this section we discuss the properties of polarized neutron matter at zero temperature 
under the presence of a strong magnetic field. Such properties are described in terms of (/ = 

10 



0, 1) Landau parameters, which are calculated from Skyrme and Gogny effective interactions. 
At each thermodynamic condition (p, T = 0,B) the corresponding induced magnetization of 
the neutron system has been obtained by minimizing the free energy within a Hartree-Fock 
calculation, as described in [l^. In order to be physically meaningful we have considered 
that the nucleon picture holds up to a maximum limiting density p = 4po (po = 0.1737/m^^). 



Similarly, more detai 
where nuclear pasta 



ed calculations are needed when considering densities below p = 0.5po 



24 



26| may be present. 



To measure the effects of the magnetic field on the Landau parameters we calculate the 
variation ratios, Rpi and Rqi as defined in Eqs.f llSI) and ( fT6l) for the dominant neutron 
component with magnetic moment aligned parallel to the magnetic field. This corresponds 
to the spin antialigned component due to the fact that /i„ < 0. 

In Fig. [1] we show the ratios Rpo, Rgq: Rpi Rgi omit the superscript) computed 
with the Skyrme SLy7 (a) and with the Gogny DIP (b) interactions for a characteristic 
low density of p = 0.5po as a function of the logarithm of the magnetic field. At this 
low density the change due to the magnetic field strength in the monopolar and dipolar 
coefficients is very mild. For Skyrme the ph interactions are mostly slightly more repulsive 
in the density and spin channel but for Gogny this behaviour is reversed. The variation of 
the dipolar terms, as computed with the Gogny interaction is larger than for the monopolar 
case. However, dipolar terms are about one order of magnitude smaller 31i so the overall 
effect is mostly unchanged with respect to the monopolar description. For magnetic field 
strengths below B ^ 10^^ G the change in the ratios is negligible due to the tiny value of the 
neutron magnetic moment, which results in a vanishing value of the induced magnetization, 
i. e., of the polarization of the neutron matter. 

In Fig. [2] we show the same ratios shown in Fig. [1] as computed with the Skyrme (a) and 
Gogny DIP (b) interactions for a high density case, p = 3.3po as a function of the logarithm 
of the magnetic field. We can see that for the Skyrme SLy7 interaction the presence of a 
ferromagnetic transition at density close to the one selected [l^, induces a larger change 
than at lower densities in the monopolar and dipolar terms as magnetic field grows. The 
attraction in the ph density excitation channel increases while in the spin channel becomes 
more repulsive as the magnetic field strength grows. For the dipolar terms this tendency 
is reversed. For the Gogny interaction the variation of the ratios is very small except for 
the Gi term. However, the large variation of Gi does not modify the tendency shown by 
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16.5 17 17.5 
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16.5 17 17.5 18 

lOQio B (G) 

FIG. 1: Variation ratio of Landau coefficients 
with the SLy7 (a) and DIP (b) parametrizations 
at density p = 0.5po and zero temperature as a 
function of the logarithm of the magnetic field 
strength. 

the monopolar terms since in this range of densities it is also an order of magnitude smaller 
than Gq. 

We now analyze some other quantities related to the Landau parameters in polarized 
neutron matter. In Fig. [3] we show the effective neutron mass for the Skyrme SLy7 (a) and 
Gogny DIP (b) interactions as a function of the logarithm of the magnetic field strength. In 
each plot the spin down (upper curve) and spin up (lower curve) polarized components at 
saturation density, po? are shown. Skyrme interaction predicts smaller effective masses and 
a larger relative variation for the spin up and down components than in the Gogny case. 
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16.5 17 17.5 

lOQio B (G) 
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16.5 17 17.5 18 

logioB(G) 

FIG. 2: Same as Fig. [1] computed at density p = 
3.3po- 

These two effects will, in turn, largely affect the level density at the Fermi surfaces of the 
spin components in magnetized neutron matter. 

In Fig. m the two spin contributions to the isothermal compressibility for neutron matter 
are shown as a function of density for a magnetic field strength B = 5 x 10^^ G and zero 
temperature for the Skyrme SLy7 (solid line) and Gogny DIP (dashed line) interactions. 
Contributions from the spin down (up) polarized particles correspond to the upper (lower) 
curves for each interaction model. In some works compressibility values are plot related to 
the Fermi gas value, in our case, due to the fact that at each density the relative populations 
of the up and down spin polarized particles change, we have chosen to plot absolute values. 
Due to the mild variation of the Landau coefficients with the magnetic field (lower than 
20%) we can see that the main effect of the magnetization in the system is the change in the 
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FIG. 3: Effective neutron mass at density po as 
a function of the logarithm of the magnetic field 
strength for the Skyrme SLy7 (a) and Gogny DIP 
(b) interactions. For each model upper (lower) 
curves correspond to spin down (up) polarized 
particles. 
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level densities. As the density increases the dominant spin down fraction of the polarized 
plasma becomes stiffer. However, there is a dramatic decrease in the global compressibility 
in the proximity of the ferromagnetic transition in the Skyrme case. This transition density 
decreases as the magnetic field grows [lO^. In fact, the results beyond this density are not 
physically meaningful but we show them for the sake of completeness in the density range 
considered in this work. In Fig. [5] contributions to the isothermal compressibility for both 
spin components in neutron matter are shown as a function of the logarithm of the magnetic 
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FIG. 4: Isothermal compressibility components obtained for the Skyrme SLy7 (solid curve) and 
Gogny DIP (dashed line) models as a function of density. Spin down (upper curve) and up (lower 
curve) polarized components are shown for each model. 

field strength at saturation density for the Skyrme SLy7 (solid line) and Gogny DIP (dashed 
line) interactions. For each model upper (lower) curves refer to spin down (up) polarized 
particles. As the magnetic field strength increases there is a splitting of the spin up and down 
component behaviour with respect to the unpolarized case (A = 0). Global compressibility 
is obtained by adding the contributions of both spin fractions and changes very mildly with 
the magnetic field. The relative variation of both components at the maximum field strength 
is of 20% (28%) for the Skyrme SLy7 (Gogny DIP) interaction at this density. 

We now consider the effect of a strong magnetic field on the static magnetic susceptibility. 
In Fig. [6] we plot the magnetic susceptibility for neutron matter as a function of density in 
units of the nuclear magnetic moment squared, /x^, for a magnetic field B = 5 x 10^^ G. The 
results obtained with Skyrme SLy7 (Gogny DIP) interaction are shown with solid (dashed) 
line. For Skyrme results we can see that the onset of a ferromagnetic transition (that for 
this specific parametrization takes place at a density around p = 3.3po and is signaled with a 
vertical dotted line on the plot) drives a divergence of the susceptibility and the subsequent 
second order phase transition. In this model the upper curve corresponds to the X~/f^n 
value and the lower curve to the x^//^n- The onset of the ferromagnetic transition is due 
to the vanishing values approached by the quantities in the denominator of Eq.( l20i) . For 
the Skyrme SLy7 interaction the meaningful regime, though, is limited for values before the 
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FIG. 5: Isothermal compressibility components obtained for the Skyrme SLy7 (solid curve) and 
Gogny DIP (dashed line) models as a function of the logarithm of the magnetic field strength at 
density pQ. Upper (lower) curves refer to spin down (up) polarized contributions. 

transition takes place. For the Gogny DIP interaction (dashed line) the splitting between the 
contributions from the two spin orientations is less than 1%, therefore, indistinguishable on 
the plot. There is no ferromagnetic transition for this parametrization and the susceptibility 
keeps finite values for the whole density range. 

400 

300 

J 200 
> 

CD 

100 

-100 
-200 

0.5 1 1.5 2 2.5 3 3.5 4 

P/Po 

FIG. 6: Magnetic susceptibility as a function of density for the Skyrme SLy7 (solid line) and Gogny 
DIP (dashed line) models. See text for details. 

In Fig. [7] we show the magnetic susceptibility in units of the nuclear magnetic moment 
squared, /x^, for neutron matter as a function of the logarithm of the magnetic field strength 
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at saturation density. We plot results obtained with the Skyrme SLy7 (Gogny DIP) with 
solid (dashed) line. Upper (lower) curves for each model refer to spin down (up) polarized 
population fractions. We see that contributions from the down (up) polarized components 
give larger (smaller) contribution to the total susceptibility that can be obtained as the 
sum of both. Notice that the global susceptibility remains almost unchanged as magnetic 
field grows. At this density the Skyrme model predicts a higher tendency of the system to 
achieve a net magnetization than in the Gogny case as B increases. No transition is near at 
this density for the Skyrme model and a smooth behaviour is observed in the full range of 
variation of the magnetic field. 



> 
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16.5 17 17.5 18 

lOQio B (G) 

FIG. 7: Magnetic susceptibility for a pure neutron system as a function of the logarithm of the 
magnetic field strength for the Skyrme SLy7 (solid line) and Gogny DIP (dashed line) at saturation 
density. For each model upper (lower) curves refer to down (up) spin polarized components. See 
text for details. 



IV. SUMMARY AND CONCLUSIONS 

In this work we have investigated in the context of the Landau Theory of normal Fermi 
Liquids, the variation of monopolar (/ = 0) and dipolar (/ = 1) Landau parameters which de- 
scribe the particle-hole interaction matrix elements for pure neutron matter in the presence 
of a strong magnetic field at zero temperature. We have used effective nuclear interactions 
such as the zero-range Skryme SLy7 and finite-range Gogny DIP and obtained their ana- 
lytical expressions valid in the case of a polarized Fermi sea. We have computed the ratios 
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of variation of these coefficients with respect to the B = case for the dominant spin com- 
ponent in the system. For the Skyrme and Gogny interactions the variation is very mild 
keeping below \R\ ^ 20% for the maximum magnetic field strength studied in this work 
B ^ 10^^ G. This is a direct consequence of the small polarization induced by such strong 
fields in the system in the meaningful range of densities considered in this work. 

We have also analyzed the effect of the presence of a strong magnetic field in some other 
static properties in the neutron plasma as deduced from the Landau Theory of Fermi Liq- 
uids. Effective neutron masses at the Fermi surface for each spin polarized component are 
calculated for both interactions from the dipolar Fi coefficients. The magnetization in the 
plasma causes a splitting in the values of the up and down spin polarized components. Spin 
down (up) polarized particles show an increase (decrease) in their effective masses with re- 
spect to the non polarized the magnetic field strength grows. This will affect the 
level densities in polarized neutron matter. In addition, Skyrme interactions predict a larger 
variation in the splitting than the Gogny forces with increasing magnetic field, B. Other 
magnitudes such as the isothermal compressibility, K, relate to the dipolar Landau coeffi- 
cients. In the presence of a strong magnetic field the compressibility is a growing function of 
density, stiffer in the Skyrme case than in the Gogny case. When a ferromagnetic transition 
is near, as it happens for the Skyrme interaction used in this work (as a representative case of 
almost all Skyrme parametrizations jl^ ) , there is a dramatic decrease of the compressibility 
around the density of the onset of the ferromagnetic instability. The magnetic susceptibility, 
X, shows a divergent behaviour at densities close to that of the onset of the phase transion 
for the Skyrme case. The contribution to the total susceptibility from the polarized popu- 
lations in the system is related to the Landau coefficients in the spin channel that approach 
values driving the system to an instability. For the Gogny case there is no divergence since 
magnetization of the system keeps very mild. For densities in the intermediate range, not 
close to the transition, spin down (up) contributions give a higher (lower) susceptibility as 
it is more (less) easy to polarize the system as the external field strengthens. Although a lot 
of work has been devoted to the study of polarized nuclear systems, additional effort should 
be made for exploring properties of nuclear asymmetric matter at finite temperature. The 
presence of magnetic fields should be carefully studied in the future with direct simulation 
techniques for the low density case. This will allow to explore the possibilities that exotic 
matter shapes in the crust of neutron stars may give as an additional source of opacity to 
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the neutrino cooling in the early stages of Supernovae cooling. 
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